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Abstract 

In the present paper we investigate the phase transition "Coulomb-confinement" 
in the Higgs model of abelian scalar dyons - particles having both, electric e and 
magnetic g, charges. It is shown that by dual symmetry this theory is equivalent to 
scalar fields with the effective squared electric charge e* 2 = e 2 + g 1 . But the Dirac 
relation distinguishes the electric and magnetic charges of dyons. The following 
phase transition couplings are obtained in the one-loop approximation: a cr u = 
e 2 rit /47r ps 0.19, a crit = g 2 crit /^ ~ 1.29 and a* crit ps 1.48. 
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1. In the present paper we consider dyons - the particles with electric e and magnetic 
g charges. As it was shown in Refs. [1,2,4-7], dyons play an essential role in physics of 
nonabelian theories (in particular, in QCD). 

The local field theory of electrically and magnetically charged particles, so called 
"QuantumElectroMagnetoDynamics" (QEMD) [8], is presented by the Zwanziger for- 
malism [9,10], (see also [11] and [12]), which considers two vector potentials A^(x) and 
B^(x), describing one physical photon with two physical degrees of freedom. Here _B M 
is the magnetic gauge potential, which is dual to the electric gauge potential A^. This 
formalism symmetrically contains non-dual and dual interactions of gauge fields with the 
corresponding currents: j^ e > ■ A and j <m - ) • B, where j e and j™ are electric and magnetic 
currents, respectively. 

Dyons are described by the field $, having charges n\e and n 2 g {ni,n 2 G Z). The 
total system of gauge fields and dyons is given by the partition function having the 
following form in Euclidean space: 

Z = J[DA}[DB}[D<&][D§ + }e- s , (1) 

where 

S = J d 4 X L(x) = S Zw (A, B) + S gf + S {matt er) ■ (2) 

The Zwanziger action Sz w (A,B) is given by: 



S Zw (A,B) = J d 4 x 



±(n.idAA]) 2 + 1 -(n.idAB]) 2 + 



(3) 



+ % -(n ■ [d A A]){n ■ [d A B}*) - Un ■ [d A B]){n ■ [d A A]*) 

where we have used the following designations: 

[A A B}^ = A^B V - A„B„ (n ■ [A A B})^ = n v {A A B)^, 

G ' piv = ~^jiv\pG\p. (4) 

In Eqs. (jSJ) and (J3| the unit vector represents the fixed direction of the Dirac string in 
the 4-space. 



The action S( matter y. 



S(matter) J d X -^(matter) (*^) 



(5) 



describes the matter fields of dyons. 

In Eq. (J2J Sgf is the gauge-fixing action. According to Ref. [11], this action is given 
by the following equation: 



Sgf = I d 4 x 



M\ 9 Ml 

— {n ■ A) + — (n • BY 



(6) 
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and has no ghosts. Eq. © contains the mass parameters Ma and M B . 

Let us consider now the Lagrangian L^ matter ^, containing the Higgs scalar dyon field 
interacting with gauge fields A^(x) and B^(x): 

L { matter)(x) = ^\V^f - [/($), (7) 

where 

V il = d il - ieA^ - igB^, (8) 

and 

[/($) = im 2 |$| 2 + ^|$| 4 (9) 

is the Higgs potential for the dyon field $. 
The complex scalar fields: 

<§ = <P + i X (10) 

contain the Higgs (0) and Goldstone (x) boson fields. 

In the Lagrangian (JZJ) the interactions, given by terms j^A^ and j^B^, contain the 
following electric and magnetic currents: 

j; = ie[^ + V^-^(V^)+], (11) 
= ig[<S> + V^-<S>(V^) + }. (12) 

In Eq. dU) we have "sea-gull" terms e 2 A 2 |$| 2 and g 2 Bf l \§\ 2 . 
Letting 

F = dAA = -(dAB)*, (13) 
F* = d A B = (d A A)*, (14) 

and the Hodge star operation on the field tensor: 

F^i/ — ^puVpaFpa-, (15) 

it is easy to see that the free Zwanziger Lagrangian (JHJ) is invariant under the following 
duality transformations: 

F <-> F*, (dAA)^(dAB), (d A A)* <-> -(d A B)*. (16) 

We also have a dual symmetry as an invariance of the total Lagrangian L(x), provided 
that electric and magnetic charges and currents transform simultaneously according to 
the following discrete symmetry: 

e->0, 9 -> -e, 

(17) 

The action (J2J), given by Eqs. @ and (jSEJ), leads to the following field equations: 

dxF^=& dxF^ = 3 J, (18) 

2 



where currents j^, j™ are given by Eqs. (fTTl [T2"j) . 



2. The Lorentz invariance is lost in the Zwanziger Lagrangian (|3*|). because it depends on 
a fixed vector n M . However, this invariance is regained for the quantized values of coupling 
constants e and g, obeying the Dirac-Schwinger-Zwanziger (DSZ) relation for dyons: 



27m i: j, 



riij G Z. 



(19) 



As it was shown in Refs. [13-18] (see also [19]), gauge theories (abelian and nonabelian) 
can be constructed in terms of loop variables. The new Zwanziger-type Lagrangian 
does not contain a fixed vector n^, but contains loop coordinates and their derivatives. 
Considering the closed loop operators: 



A(C) = Pexp 



(20) 



which measures magnetic flux through C and creates electric flux along C, and: 



B(C) = Pexp 



ig & A^dr] 11 



c 



(21) 



which measures electric flux through C and creates magnetic flux along C, we have the 
t'Hooft commutation relation [20]: 



A{C)B(C) = B(C)A(C)exp(2mn), 



(22) 



where n is the number of times C winds around C . In Eqs. fl20j) and (|21j) we have 
coordinates £ M and of the loops C and C in the 4-dimensional space. 

Eq. (}2*2*|) gives the Dirac relation - the charge quantization condition for a single 
dyon, interacting with fields A^ and B^. 



eg = 27m, n G Z. 

For elementary particles we have: 

1 

eg = 2tt, or aa = — , 
where a and a are the electric and magnetic fine structure constants: 



(23) 



(24) 



a = — and 

47T 



a 



An' 



(25) 



3. The effective potential in the Higgs model of scalar electrodynamics was calculated in 
the one-loop approximation for the first time by authors of Ref. [21]. A general method 
of calculation of the effective potential is given in the review [22]. Using this method, 
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we can construct the effective potential for theory described by the partition function (JTJ) 
with the action S, containing the Zwanziger action (J3J), gauge fixing action @ and the 
action (JTJ) for dyon matter fields. 
Let us consider now the shift: 



(26) 



with <3>b as a background field, and calculate the following expression for the partition 
function in the one-loop approximation: 



5S($) 



+ h.c. 



<i>=<e> s 



[DA][DB] [D$][D<$>+] 

•exp -J d 4 x 

= exp{-F($ 6 ,e 2 ,# 2 ,m 2 ,A)}. 
Using the representation we obtain the effective potential: 

V eff = F(0 6 ,e 2 ,# 2 ,m 2 , A), 
given by the function F of Eq. (|2*7J) for the constant background field: 

&b = 4>b = const. 



(27) 



(2f 



(29) 



1 m 2 



The effective potential (|2*%|) has several minima. Their position depends on e 2 , g 
and A. If the first local minimum occurs at 0& = 0, it corresponds to the so-called 
"symmetrical phase" , which is the "Coulomb" phase in our description. 

We are interested in the phase transition from the Coulomb phase "06 = 0" to the 
confinement phase "06 = <po ^ 0" . In this case the one-loop effective potential for dyons 
is similar to the expression of the effective potential, calculated by authors of Ref. [21] for 
scalar electrodynamics and extended to the massive theory in Ref. [23]. 

As it was shown in Ref. [21], effective potential can be improved by the consideration 
of the renormalization group equation (see also the review [22]). 



4. The renormalization group (RG) describes an independence of a theory and its cou- 
plings on an arbitrary scale parameter M. We are interested in RG applied to the effective 
potential. The renormalization group equation (RGE) for the effective potential means 
that the potential cannot depend on a change in the arbitrary renormalization scale pa- 
rameter M: 

(30) 



dV ef f 



dM 

The effects of changing it are absorbed into changes in the coupling constants, masses and 
fields, giving so-called running quantities. Knowing the dependence on M 2 is equivalent to 
knowing the dependence on <p 2 = 2 . This dependence is given by RGE. Considering the 
RGE improvement of the potential, we follow the approach by Coleman and Weinberg [21] 
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and its extension to the massive theory [23]. Here we have the difference between the scalar 
electrodynamics [21] and scalar QEMD. 

RGE for the improved one-loop effective potential can be given in QEMD by the 
following expression: 

M ii +> 3 4 + ei3 4e + + - 7 4) v *" (,p2) = °- (3l) 

where the function 7 is the anomalous dimension: 



ilW—m- (32) 



RGE (|3*T|) leads to a new improved effective potential [21]: 



K//(0 2 ) = \m 2 ren {t)G\t)<f + \\ ren {t)G\t)<p\ (33) 



where 



G(t) = exp 



1 " 



dt'j(g ren (t'), X r en(t')) 







2 

Eq. (JBTj) reproduces also a set of ordinary differential equations 

d\ rP/n 



(34) 



(3\(g r en(t), Ken(t)), (35) 



= raren(*)A™ 2 )(SVen(£), Aren(*)), (36) 
= "^T^ = ^en(*), W*)) - P 9 (9ren(t), A ren (t)) = (37) 

where t = ln(0 2 /M 2 ), and the subscript "ren" means the "renormalized" quantity. 

The last equation (J37)l is obtained with the help of the Dirac relation (|24|) for minimal 
charges. Indeed, in Eq. (|3*T|): 



9 ^ d d Q de d / fl / 2vr\\ 9 

ae ae dg oe \ \ g 2 J J oe 



^e-^T^rz = ^ total) T^rz, (38) 



_^ _ p(total) ^_ 

d In e (9 In e 

where 

We can determine both beta functions for X ren and rr? ren by considering a change in M in 
the conventional non-improved one-loop potential. 
Let us write now the one-loop potential: 

V eff = V + V u (40) 
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where 

Vo = ^ + \<t>\ (41) 

and (see Refs. [21] and [22]): 



64vr 2 



^2 ^2 i q\a2 



3(e 2 + /) 2 4 ln^- + (m 2 + 3A0Tln 



2l m 2 + 3A0 2 



M 2 r y M 2 



+ (m 2 + A0 2 ) ln^^ 2m 4 In — 



(42) 



We can plug V e ff, given by Eqs. ([4UH42jL into RGE (j3~Tj) and obtain the following equation 
(see [22]): 

Equating 2 and 4 coefficients, we obtain: 

Pi \2 o( 2 , 2 \2 

/3 Q~ \ i ren i V ren 1 Jren) I A A\ 

= 2lKen + ^ + 16^ ' (44) 

/V)=7 + ^- (45) 

The result for 7 is given in Ref. [21] for scalar field with electric charge e, but it is easy 
to rewrite this 7-expression for dyons with renormalized charges e ren and g ren : 

qf p 2 , 2 \ 

V ren 1 iJren) (46) 

167T 2 



Finally we have: 



^ ~~ 167T 2 [^^ 6ren ^ ren ^ ^^ren 6^ren(& r en + 5Ven)] ) (47) 



ren ren 



&l ,, ; 2 [4A 2 en -3(e 2 en + ^ 2 e J]. (48) 
Now the aim is to calculate the beta-function p( total ) in Eq. (|37j) . 

5. Many years ago it was calculated (see Refs. [24-28]) that Gell-Mann-Low equation [29]: 

din a „. . . . 

— =/?(«) (49) 

has the following /5-function for the electric charge in the scalar electrodynamics: 

« o) = is( 1+3 s + -)- (50) 
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According to the previous item 4, now we have: 



VE£L = -*E* = f, (a) - m = fP~> (a) . (51) 

These RGEs are a consequence of the Dirac relation and the dual symmetry, consid- 
ered in the item 1. 

If both a and a are sufficiently small, then beta-functions in Eq. ()51j) are described 
by the contributions of the electrically and magnetically charged dyon loops simultane- 
ously Their analytical expressions are given by (|5()jl . and in the two-loop approximation 
we have the following equations (|51|) for scalar dyons (see also [12]): 

dlna dln& a-5/. . 3 « + a | \ = fitotai),^ ^ 



dt dt 12vr V 4?r 

According to Eq. the two-loop contribution is not more than 30% if both a and a 
obey the following requirement (see Ref. [12]): 

0.25<a,3<l. (53) 

The lattice simulations of compact QED give the behaviour of the effective fine structure 
constant in the vicinity of the phase transition point [30-34] . The following critical values 
of the fine structure constants a and a were obtained in Ref. [30]: 

affit « 0.20 ± 0.015, a l £ it « 1.25 ± 0.10. (54) 

Eq. (jHljl demonstrates that a cr u and a cr it, obtained in the compact lattice QED [30-34], 
almost coincide with the borders of the requirement (|53|h given by the perturbation 
theory for /^-function [12]. Assuming that in the vicinity of the phase transition point 
the coupling constant g ren may be described by the one-loop approximation, we obtain 
from Eq. (l52f the following RGE for g 2 



) ren' 



dt 48tt 2 
Using the Dirac relation ()24j) . we have: 

do 2 o 4 1 

tiren jren 

dt ~ 48^ ~ ^2' 



d\n(o 2 ) a 2 - e 2 

1 ■> nu ■ ■inn fen (55) 



(56) 



Note that the second term of Eq. (|56|l describes the influence of the electric charge on the 
behaviour of the magnetic one. 



6. In this part of our paper we investigate the phase transition from the "Coulomb" phase 
(4>b = 0) to the confinement one (</>& = O 7^ 0), following the methods of Refs. [35-37]. 
This means that the effective potential (J3*3*j) of the Higgs scalar dyons has the first and 
the second minima appearing at <p = and cf) = <po, respectively. They are shown in 
Fig. dby the solid curve "1". These minima of V e ff{4> 2 ) correspond to the different vacua 
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arising in this model. The conditions for the existence of degenerate vacua are given by 
the following equations: 

Veff(P) = V ef f(<ft) = 0, (57) 



dV f 



eff 



dV f 



eff 



0=0 



<P=<f>0 



with inequalities 



d 2 V e 



eff 



> o, 



d 2 V f 



eff 



<f>=0 



> 0, 



4>=<f>o 



or considering <p 2 as a variable, we can write: 

_ dV eff 



V "fMD 



4>=4>o 



d 2 V f 



eff 



d{<p 2 



> 0. 



(58) 
(59) 

(60) 
(61) 



From now we omit (for simplicity) the subscript "ren" , using the following designations: 



2 

ren' 



\ 2 = \ 

9i "j = &reni 9reni ^reni m> 

The first equation (JSTj) applied to Eq. (J3~3~j) gives: 

m 2 = --A^o^o), where t () = h, ( j± 



(62) 



The calculation of the first derivative of K//(0 2 ) leads to the following expression: 



K//(^ 2 ) 



V e ff(4> 2 ) A . „dlnG\ ldm 2 



1 + 2- 



dt 



2 rft 



G 2 (t) 



1 A , . d\ , dlnG\ | , 
+4 A W + * +2A — l G( ^ 



(63) 



From Eq. an d ()46|) we have 



dlnG 1 3(e 2 + c/ 2 ) 



dt 



"2 7 



32tt 2 



(64) 



Using Eqs. (}4Tj) . (JUJ, (|52*|) and (|54*j) . it is easy to find the joint solution of equations 
VeffM = K//(0o) = 0: 



or 



VZffM = t + A + ft - 7 A) G\t ) 



(3 X + A (1 - 7 - /3 {m2) ) = 0. 



^o = 0, 



(65) 
(66) 
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Eq. (job)) gives the phase transition border, valid in any loop approximation. Putting into 
Eq. (pjj) the function /3\, /3( m 2) and 7, given by Eqs. (j44H46|) and (foljj) in the one-loop 
approximation, we obtain the following equation for the phase transition border: 



V 



effWo) 



or 



16?r 2 
3 



0. 



Using the Dirac relation 



4 

9crit 



167T 2 

and Eq 

4 



(2tt, 
+ + 8tt 2 



9crit 



, we have: 
-2A 



57T 



+ A 



(67) 
(68) 

(69) 



The curve ()69|) is shown on the phase diagram (A; g 2 ) of Fig. [21 by the curve "1", which 
describes a border between the "Coulomb" phase with V e ff > and the confinement ones, 
having V e J f n < 0. 



7. The next step is the calculation of the second derivative of the effective potential: 
„„ , l2s V' e ffU 2 ) ( 1 , ld 2 m 2 dm 2 dlnG 2 d 2 \nG 

, /dlnGA 2 \ G 2 /IdA ld 2 A c/ArflnG dhxG 
+2m 2 — — —+ _— + -— + 2—— ^ + 2A- 



tft J <p 2 V 2 rf* 4 rft 2 rft tft tft 
, d 2 lnG ^ /dlnG\ 2 \ . 

+ A — +4A (— ) J G w- < 70 > 

Let us consider now the case, when this second derivative changes its sign giving a max- 
imum of Veff, instead of the minimum at <fi 2 = 0q. Such a possibility is shown in Fig. ^ 
by the curve "2" . Now the two additional minima at <fi 2 = 2 and 2 = (f)\ appear in our 
theory. They correspond to the two different confinement phases related with the con- 
finement of the electrically charged particles. If these two minima are degenerate, then 
we have the following requirements: 

VefM) = Veff^j) < (71) 

and 

V'efM) = V'efM) = 0, (72) 

which describe the border between the confinement phases "conf. 1" and "conf. 2", pre- 
sented in Fig. |21 This border is shown by the curve "3" at the phase diagram (A; g 2 ) of 
Fig. El The curve "3" meets the curve "1" at the triple point A. It is obvious that, accord- 
ing to the illustration of Fig. ^ this triple point A is given by the following requirements: 

VefM) = Kfti&) = VefMl) = o- (73) 
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In contrast to the requirements: 

VefM) = v; ff (<ft) = o, (74) 



VefM) = K/ f (0o) = 0. (75) 



describing the curve "1", we are going now to consider the joint solution of the following 
equations: 

It is possible to obtain this solution using Eqs. flZOJ), (JE21), (03), AH and |fl2$: 

F(\,g 2 )-7i 2 \ + 27i 2 (e 2 + g 2 ) = 0, (76) 

where 

jF(A,g 2 ) = 5(e 2 + ^) 3 + (247r 2 + 12A)(e 2 + ^) 2 -9A 2 (e 2 + ( ? 2 ) 

+36A 3 + 80vr 2 A 2 + 64vr 4 A. (77) 

According to the Dirac relation, we have e 2 = (2ir/g) 2 and obtain the curve "3" of Fig. 
which represents the solution of Eq. ()76|) (together with 1)770). which is equivalent to 
Eqs. 013). 

The intersection of the curve "1" (or "2") by the curve "3" gives the value of the 
triple point A (or B), which corresponds to the critical values a cr u and a cr u in the Higgs 
model of scalar dyons. 

The numerical calculations demonstrate that the triple point A exists in the very 
neighbourhood of the maximum of the phase transition curve "1" of Fig.0 and its position 
is given by the following values of A and g 2 : 

\ A « -13.44 and g\ « 16.16. (78) 

The triple point B has its position near the minimum of the curve "2" and gives the 
following values of A and g 2 : 

\ B « -13.44 and g\ « 2.44. (79) 

Here we see that the ratio: 

e 2 A 47r 2 . , 

^ = - r «0.15 (80) 
9a 9a 

describes 15% contribution from e\ to the value g\. This result explains the difference 
between the dyon model and scalar electrodynamics (see the review [37]). 

The negative values of \a,b in the triple points are not dangerous, because they are 
obtained for the renormalized quantity, and point out that <f) Q M cuto ff, what is quite 
possible. 

It seems to us, that the triple point A (or B) does not coincide with the maximum 
(minimum) of the curve "1" ("2") due to our one-loop approximation. We expect that 
they will coincide in the true theory. 

Values (|7HD and JSJ) for g\ 

b correspond to the following critical fine structure 

constants: 

a crit ^0.19, and a crtt » 1.29, (81) 
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which are in agreement with the result of lattice investigations for QED. 

In the present theory of scalar dyons, the expressions (JHEJ), (EH) and (J77J) depend on 
the quantity e* 2 = e 2 + g 2 . This is a consequence of the dual symmetry (|16|) and (|T7j) . 
By dual transformations, we can reduce dyon fields to the electric fields with the effective 
electric charge e*. From Eq. (|8~T]) the critical value of a* = e* 2 /4n is: 

<nt « 1-48. (82) 

But the Dirac relation ((23), obtained in the item 2, distinguishes the electric and magnetic 
charges of dyons, when it interacts with the fields and B^. This is significant for the 
phenomenon of the string formation in QCD, considered in Ref. [7]. 
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Fig. 1: The effective potential V e ff the curve "1" corresponds to the "Coulomb- 
confinement" phase transition; curve "2" describes the existence of two minima corre- 
sponding to the confinement phases. 
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Fig. 2: The phase diagram (A; g 2 ) for the Higgs model of scalar abelian dyons. Curves "1" 
and "2" separate "Coulomb" phase from confinement phase. Figure shows the existence 
of triple points A(X A = -13.44; g\ = 16.16) and B(\ B = -13.44; g% = 2.44). These 
triple points are boundary points of three phase transitions: the "Coulomb" phase and 
two confinement phases: "conf. 1" and "conf. 2", which are separated by the curve "3". 
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